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first two peaks grow with temperature. The fact that these two peaks also correspond 
to electric quadrupole transitions is explained by a recently proposed theory. 
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1. Introduction 

Pre-edge peaks often arise at the K-edge of transition metal elements. This pre-edge 
structure is sensitive to the metal valence, to the symmetry of its surroundings and to 
the atomic species of the neighbours (see [1] for a recent review). As a consequence, the 
measurement and analysis of the pre-edge peaks are widely used in earth sciences [2], 
biology [3] , chemistry [I] and physics [51 [6] . 

Because of their practical importance, pre-edge features have to be well understood 
and they were the object of detailed theoretical works using various approaches: 
multiplets [3 IB], Bethe-Salpeter equation [9], multiple-scattering [in], pseudopotential 
theory [51 [11] . Vedrinskii and his group were particularly active to extract information 
from the pre-edge structure [121 131 EH] • 

In section 2, we make a short review of the literature to show that the temperature 
dependence of pre-edge peaks is not a rare property of x-ray absorption spectra. 
However, this dependence is usually attributed to static off-center displacements or to 
phase transitions. Therefore, our preliminary investigation [15] showing a temperature 
dependence of the pre-edge peaks at the titanium K-edge in Ti0 2 (rutile) came 
out as a surprise because the pre-edge variation was observed in a temperature 
range where no phase transition occurs and where many high-precision structural 
studies [El H71 HB1 EIHl ED] indicate that no off-center atomic displacement takes place. 
Soft modes have indeed be reported [21 J but the calculated phonon spectrum shows 
excellent agreement with experiment and no imaginary mode is present j22j [231 EH 125] 
when the proper functionals are used [26J. 

Thus, we carried out detailed experiments to confirm and analyze this temperature 
dependence. The results of these experiments are presented in section 3. Section 4 
describes why such a temperature dependence is a priori surprising and sketches a 
theoretical interpration that enables us to understand why the temperature dependence 
is restricted to the first two peaks and why no energy shift is observed. A conclusion 
summarizes our results and provides possible extensions of this work. 

2. A short review 

In this section, we present a short and non-exhaustive review of the temperature 
dependence of pre-edge peaks. 

As far as we know, such a temperature dependence was first observed by Durmeyer 
and coll. [15] at the K-edge of titanium in Ti0 2 (rutile), Li 4/ / 3 Ti 5/ / 3 04 and LiTi 2 4 . 
It was subsequently measured at the titanium K-edge of a several perovskite crystals: 
PbTi0 3 [23 EB1 CE21 ES ED], SrTi0 3 [HI ED], BaTi0 3 [30], CaTi0 3 [30]. 

A similar temperature dependence was observed at other edges in perovskite 
crystals: at the niobium K-edge of KNb0 3 [321 [33], NaNb0 3 [31], PbIni /2 Nbi /2 3 [35], 
at the zirconium K-edge of PbZr0 3 [14|, PbZr 515 Ti .485O 3 [H], BaZr0 3 [14J, at the 
K-edge of Mn in Lai_ x Ca x Mn0 3 [36], [37] and at the K-edge of Fe in in La .8Sr . 2 FeO 3 
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and Lao.7Sr .2Fe0 3 [38]. In most cases, the temperature effect was intepreted in terms 
of a phase transition or of a static off-center atomic displacement due to the presence 
of very soft modes in the crystal. 

However, the effect is not restricted to the perovskite structure. Apart from the 
results of Durmeyer and coll. [15J, it was observed at the K-edge of titanium in Ti02 
and Mg 2 Ti04 [30], at the L-edges of La in Sr-doped La 2 Cu0 4 [39J and at the K-edge 
of V in V0 2 [ID]. A temperature dependence of XANES spectra was also observed at 
the K-edge of oxygen in water [H] and in doped LaMnOs [^2l |4"3] . Finally, the Mahan- 
Nozieres-Dominics singularity can also give rise to a temperature dependence of the 
x-ray absorption spectra of metals (see Ref. |Sj for a review). 

We come now to our experimental temperature dependence at the K-edge of 
titanium in rutile. 

3. Experiment 

The x-ray absorption experiments were performed at the Dll (energy dispersive) and 
at the EXAFSII stations of the DCI storage ring of the Laboratoire pour l'Utilisation 
du Rayonnement Synchrotron in Orsay (France). 

A rutile single- crystal plate (9 mm x 4 mm x 50 /jm) was measured at the Dll 
station in the transmission mode. The crystal plate was placed inside a liquid-helium 
cryostat operating between 4.2 K and 300 K. Measurements were carried out for two 
orientations, with the (110) face of the crystal perpendicular to the x-ray beam and 
the c-axis either parallel or perpendicular to the linear polarization vector of the beam. 
The polychromator consisted in a curved Si(lll) crystal focusing the X-ray beam at 
the center of the cryostat sample holder. Higher harmonics were rejected by a Si0 2 
plane mirror. The x-ray intensity was measured by a photodiode array detector. Each 
spectrum was obtained as a result of four measurements: Iq (without sample and with 
the beam), /obiack (without sample and without beam), / (with sample and with the 
beam), /black (with sample and without beam). The absorption spectrum was then 
obtained from the formula a = log(/o — /obiaek) — log(/ — /black)- The x-ray energy 
corresponding to each detector pixel was determined by comparing the spectra with a 
spectrum measured on a two-crystal monochromator beamline. The energy resolution 
was typically 0.8 eV. 

Our preliminary study [15] showed us that the pre-edge structure could exhibit a 
low signal to noise (S/N) ratio when the crystal thickness was optimized for the edge- 
jump. Therefore, we optimized the crystal thickness for the pre-edge structure and cut 
an approximately 50 fim thick crystal plate. As a consequence, we obtained excellent 
spectra in the pre-edge region but the XANES spectra after the edge had a rather 
low S/N ratio and, for each polarization direction, we normalized the spectrum at the 
inflection point of the absorption edge instead of at the edge jump. 

To check the validity of this procedure, we carried out additional experiments 
at the EXAFSII station. The experimental equipment consisted of a two-crystal 



Temperature dependence of pre- edge XAS peaks 




4965 4970 
Energy (eV) 



Figure 1. Temperature variation of the pre-edge peaks at the K-edge of titanium in 
rutilc, with the x-ray polarization vector parallel and perpendicular to the c-axis of 
the crystal. 



Si(311) monochromator, an ionisation chamber to measure the incident beam and an 
electron-yield detector. We measured a bulk rutile single crystal with the (110) face 
perpendicular to the x-ray beam and with the c-axis of the crystal either parallel or 
perpendicular to the x-ray polarization vector. The S/N ratio of the pre-edge region was 
comparable to that of the Ti K-edge spectra of rutile measured on the same beamline in 
similar conditions [431 14~6] . The experimental spectra were normalized by the standard 
procedure and, as in our previous work [15J, the temperature dependence was found 
to be negligible except in the pre-edge region. Moreover, the observed spectra and 
temperature dependence agreed well with the transmission experiments of the Dll 
station. In the present paper we show only the results of the transmission experiments 
because of their better S/N ratio. 



4. Experimental results 

Figure [H shows the pre-edge features of rutile recorded at different temperatures with 
the polarization vector perpendicular and parallel to the crystal c-axis. A decrease of 
the first two peaks Ai and A2 is observed at low temperature, whereas the third peak 
A3 does not show any significant variation. It is important to notice that temperature 
induces a change in the intensity but not in the energy position of the peaks. 
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Figure 2. Relative intensity of peak as a function of temperature for e\\c, 
normalized to 1 at 8K. Dots: experimental results; thick solid line: fit to the function 
a(l + c~ e/T ) + (1 - a); thin solid line: fit to the function acosh(6>/2T) + (1 - a). 

The physical origin of the pre-edge peaks of titanium in rutile is well known |47j . 
Peaks Ai and A 2 correspond to electric quadrupole transitions towards 3d states of 
titanium with t 2g and e g symmetry, respectively. Therefore, the peaks that vary with 
temperature are also the peaks corresponding to quadrupole transitions. 

In the optical range, the effect of temperature is usually described by a simple 
model developed by Holmes and McClure [151 HHJ ED], in which the intensity of the 
vibronic peak varies as 1 + e~ 9 ^ T , where 9 is the energy of the first vibrational level. 
Figure [2] shows the variation of the A 2 peak with temperature, fitted to the function 
a(l + e- 9/T ) + (1 - a), where a(l + e S / T ) represents the fraction of the A 2 peak that is 
purely vibrational and 1 — a the fraction that is due to electric quadrupole transitions 
(and to the possible tail of the electric dipole peak A 3 ). 

The result of the fit is a = 0.21 and 9 = 168 K ±10 K. Note that the value of 9 
compares favorably to the energy of the first odd vibrational level at the V point obtained 
by ah initio calculations (168 K [51], 150 K [26], 169 K [25] or 181 K |23j) or by neutron 
scattering 163 K [52]. However, the simplicity of the Holmes and McClure model implies 
that the quality of this agreement is probably fortuitous. Indeed, an alternative single- 
mode model is sometimes used [53J [511 ESI [56] , for which the temperature dependence 
is coth#/2T. For this second model the fit gives 9 = 58 K ±5 K and a = 0.014. 

5. Interpretation 

It remains to understand why only the first two peaks vary with temperature while the 
rest of the XANES spectrum remains constant. We first describe the arguments that 
are usually given to explain the absence of temperature dependence of XANES spectra. 
Then, we show why, in some circumstances, this independence can be broken. 
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5.1. The temperature independence of XANES spectra 

There are many reasons to believe that the pre-edge features of x-ray absorption spectra 
do not depend on temperature in the absence of structural transition. The first reason 
comes from the temperature dependence of the EXAFS part of x-ray absorption spectra, 
which is represented by a Debye- Waller factor e~ 2k2<T in the EXAFS formula. The Debye- 
Waller factor accurately describes the temperature dependence of XAS in crystals, 
although it has to be supplemented with higher order cumulants in disordered materials. 
Moreover, it is well understood because it can be calculated ab initio with a good 
agreement with experiment [57J [58j EH [601 [61]. If we use this factor to describes the 
temperature dependence near the edge, we must take an energy very close to the Fermi 
energy, so that k is very small and the factor is close to unity. 

Of course, near the edge, the effect of temperature is not supposed to be described 
by a Debye- Waller factor and we must use a more sophisticated approach. Natoli's 
rule [62] gives good results near the edge. However, it describes an energy shift through 
the equation kR =Constant, and we do not observe any energy shift. More elaborate 
theoretical analyses were carried out. Brouder and Goulon [63], [63] used Lie group theory 
to describe the influence of a displacement on the multiple-scattering operator. However, 
the temperature evoluation is expected to be small near the edge, essentially because 
of Natoli's rule. Poiarkova and Rehr [57] extended the Debye- Waller factor to multiple- 
scattering paths. Their formalism is not really valid in the pre-edge region, but if we 
try to extend it we find a very small temperature dependence because of the presence of 
the k 2 factor in the exponent. Fujikawa [65l|66j used Schwinger's technique to calculate 
the effect of the Franck-Condon factors on XAFS. He concluded that this effect was not 
important. In a later work [67], he investigated the effect of temperature through the 
Keldysh approach to nonequilibrium systems. He found that thermal vibrations could 
be represented by a convolution with the phonon spectral function. His result is valid in 
the pre-edge region but leads to a small temperature effect. Moreover, the convolution 
should give rise to a broadening of the pre-edge peaks with temperature. Again, this is 
not compatible with our experimental results. A further elaboration of his approach [68] 
lead to similar results. 

We can try to take vibrations into account by coming back to the Born- 
Oppenheimer approximation and writing the initial and final wavefunctions as a product 
of a vibrational and an electronic function. However, this approach looks like a dead 
end if we consider the work by Mader and Baroni [69J who showed that, at the K-edge of 
carbon, the vibrations in the final state are strongly anharmonic and are deeply affected 
by the presence of the core-hole. Therefore, we are not allowed to consider the vibrations 
as similar in the initial and final states and we cannot use the harmonic approximation. 

Ankudinov and Rehr [70] brought some hope by showing that the S K-edge 
spectrum of SF 6 is closer to experiment when the atomic positions are slightly shifted 
with respect to the equilibrium position. But, as can be seen in their figure, atomic 
displacements shift the position of the main lines and this shift is not experimentally 
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observed. 

5.2. The temperature dependence of XANES spectra 

We recently proposed a model that enables us to understand the observed temperature 
dependence [TTj . Although a detailed account of this model would be beyond the scope 
of the present paper, we can give a physical description of the underlying physics. 

We start from the Born-Oppenheimer approximation where the wavefunctions of 
the electron+nuclei system is the product of a vibrational function by a solution of the 
Schrodinger equation for clamped nuclei. The energy of these wavefunctions do not 
depend on the position of the nuclei (as the eigenvalues of the Schrodinger equation 
for an electron in a potential do not depend on position). The transitions are made 
between these wavefunctions. If we assume that the vibrational energies are small with 
respect to experimental resolution, we can sum over the final state vibrational functions 
and we obtain an average over the vibrational function of the initial state of transitions 
for which the transition energy do not depend on the atomic positions. This explains 
why the peak positions do not move while they move if we calculate the spectrum of a 
distorted structure. 

The second step of the model consists in making a different approximation for the 
initial and the final states. The initial state is taken to be the core state centered at 
the position specified by the vibrational wave function. For the final state, we make the 
crude Born-Oppenheimer approximation, where the electronic wavefunction is assumed 
independent of the position of the absorbing atom. Then, the cross-section boils down 
to an average of the x-ray absorption spectra for a shifted core wavefunction (with 
fixed energies). What happens next can be sketched by an oversimplified model of the 
shifted core wavefunction. We assume that the displacement R is small compared to 
the electronic variable r and we obtain, to first order in R and for a spherical core state 
<fio (r), the shifted function 



When multiplied by e • r, the additional term gives us a factor e • rR • r that can be 
transformed into the sum of a monopole term proportional to (e • R)r 2 and a quadrupole 
term. The monopole term gives rise to transitions towards s states, the quadrupole 
term to transition towards d states. The transition towards s states are observed at the 
aluminium or silicon K-edge [71] . the transitions towards d states are observed at the 
K-edge of transition metals because of the presence of a strong density of d states. This 
explains why the temperature variation occurs at the position of the quadrupole peaks. 
Finally, the fact the temperature-dependent pre-edge peaks grow with temperature is 
due to the corresponding increase in thermal vibration amplitudes. 

Of course, eq. (PQ) is not sufficient because the integration over r includes also a 
region where r < R. The full theory [71] is more complex but the physical idea is the 
same. 




r R 



(1) 



r 
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6. Conclusion 

In this paper, we have presented the temperature-dependence of pre-edge features at 
the K-edge of titanium in rutile. This temperature dependence is not due to a phase 
transition or to a static distortion of the titanium site. 

The temperature dependence changes only the intensities of the peaks and not their 
positions. Moreover, the peaks that vary with temperature are the electric quadrupole 
peaks of the spectrum. An explanation of this behavior was given in terms of the 
dynamic displacement of the absorbing atom. 

Two conditions turn out to be crucial to observe temperature-dependent pre-edge 
peaks at the K-edge : (i) a large density of d states below the p states, so that 
the transitions to final d states are significant and visible; (ii) the existence of low- 
energy vibrational modes, so that the temperature effect can be observed at reasonable 
temperatures. Both of these conditions are satisfied in rutile and in perovskites 
containing transition metals. In that case, the temperature dependence provides 
information on the local vibrations around the absorbing atom. This can be particularly 
useful to investigate the vibrations of transition metal impurities. 
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